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Abstract 



QQ ' At the classical level anyons with harmonic confinement are known to 

O . exhibit two important properties namely partial separability and pseudo- 

Qs ' integrability. These stem from the fact that this system is locally identi- 

"7^ ' cal to isotropic oscillator system but differs in the global topology of the 

■^ ' phase space. We clarify the meaning of pseudo-integrability and show that 
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^ ' it amounts to a definite reduction of the symmetry group. We elaborate 

on the role of the fundamental group of the phase space and pseudo- 

r> I intrgrability in the context of periodic orbit theory and obtain evidence of 



non-exactly known eigenvalues from the semiclassical trace formula. We 
also discuss an ambiguity regarding the 'half period' trajectories suggested 
by classical modeling and exhibited by the exactly known propagator for 
two anyons. 
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I. INTRODUCTION 

Anyons as a research area is now over twenty years old |I[]. These systems have 
emerged as being interesting in their own right from mathematical physics point of view 
both at the classical and the quantum level. These systems constitute an example of in- 
equivalent quantizations due to non trivial fundamental group of the configuration space 
0. These two dimensional systems have a configuration space whose fundamental group 
is a braid group [§]. Its one dimensional unitary representations give the kinematic clas- 
sification labelled by the so called "statistical" parameter a. Being a system of finitely 
many degrees of freedom these are simpler to analyze for the dynamical consequences 
such as the spectrum of the Hamiltonian. 

The quantum mechanical spectrum (with harmonic potential added) shows two dis- 
tinct qualitative features: a) eigenvalues values which depend linearly on the statistical 
parameter, a, all of which are exactly known and b) eigenvalues which depend non- 
linearly on a and none of which is exactly known. These features were traced to the 
properties of partial separability and pseudo-integrabihty manifested at the classical 
level @. 

Partial separability of the Hamiltonian was exhibited explicitly in terms of two col- 
lective degrees of freedom and the remaining "relative" degrees of freedom. The to- 
tal Hamiltonian, after removing the center of mass degrees, was obtained as a sum of 
Hi{collective) + H2{collective, relative) . The commutator of Hi and H2 was shown to 
be proportional to if 2- This implies that the subspace of eigenstates states of the full H 
on which H2 vanishes are exact eigenstates states of Hi and these give all the exactly 



known eigenvalues values and eigenstates states. 

Pseudo-integrability, a concept first introduced by Ricliens and Berry 0, was de- 
scribed by exhibiting 2N constants of motion in involution which fail to lead to inte- 
grability via usual action-angle coordinates. This was conjectured to be the reason for 
the level repulsion seen numerically in the non-linearly interpolating spectrum |^. The 
precise meaning of pseudo-integrability in this context, however, was not elaborated. In 
particular, how this classical feature translates into there being only two good quantum 
numbers, was not analysed. 

There is another approach one could follow to get a handle on the many anyon spec- 
trum namely the stationary phase approximation (SPA) to the propagator G{E + ie) 
developed using a suitable path integral representation 0. In this approach the propa- 
gator is typically obtained as a sum over (families of) periodic trajectories in the classical 
phase space. This is entirely given in terms of classical quantities. One therefore expects 
to see directly the effects of non-trivial fundamental group of the phase space, pseudo- 
integrability and partial separability in a semiclassical framework. Since at the classical 
level, anyons with oscillator confinement, are locally identical to the oscillator system 
(but differ in the global topology of the phase space), the periodic orbits are known. 
Thus application of the periodic orbit theory (POT) to anyons, may be expected to be 
tractable. 

There is one point to be noted at the outset. At the quantum level, the statistical 
parameter, aq, enters via the stipulation of multivaluedness of the wave functions and 
is dimensionless. At the classical level its counterpart, ac-, enters as the coefficient of a 



total derivative term in the Lagrangian and has the dimensions of h. These two must 
be related as ac = hag. In the SPA computation, ac is held fixed with h going to zero. 
In the leading approximation, the spectrum will depend linearly on a^ and thus also on 
aq. Alternatively, for comparison with the quantum spectrum to leading order, one will 
want to keep a^ fixed implying that ac to go to zero with h. Viewed either way, one 
may not expect to see non-linear dependence on the statistical parameter at the level 
of the leading approximation. Nonetheless, if the non-linearly interpolating eigenvalues 
also have a linear piece, one should be able to see it in the leading SPA. 

In this work, we demonstrate that the classical analysis is adequate for exhibiting 
quantum symmetries. In this process we sharpen and clarify the meaning of pseudo- 
integrability and role of fundamental group of the phase space. 

We also sketch application of the periodic orbit theory to many anyons. In the pro- 
cess we show the presence of eigenvalues which are potentially non-linearly interpolating. 
Such eigenvalues have been seen in the numerical spectrum for three and four anyons 
1^ and have been conjectured to be present for general A^ Q]. We obtain these system- 
atically from POT. 

Further, both the classical modeling and the exact propagator for two anyons, indi- 
cate an ambiguity regarding possibility of closed orbits with half the basic period and 
its inclusion in the trace formula. We discuss this ambiguity in some details. 

The paper is organized as follows. 



In section II we specify the classical model. The system is shown to be identical to 
just the isotropic oscillator locally but not at the global topological level. The set of 
classical trajectories is detailed. 

In section III we classify the classical trajectories using symmetry transformations 
and explain how the non-trivial global topology reduces the symmetry group. 

In section IV we sketch the POT application to anyons and demonstrate existence 
of eigenvalues other than the exactly known ones. 

Section V contains a discussion of the issue of half period trajectories. We show that 
it is possible to regularise the point-like statistical flux and deduce the existence of half 
period trajectories. We look for their evidence in the exact propagator for two anyons 
and point out their ambiguous status. 

Section VI contains general remarks on the role of the fundamental group, possible 
reduction of symmetries due to non-trivial topology of the phase space. A summary of 
results and conclusions is presented. 

Appendix A contains details of 0Sp(4N, R) symmetry. 

Appendix B contains discussion of regularised classical dynamics of anyons and the 
presence of half period trajectories. 



II. CLASSICAL MODELING AND DYNAMICAL TRAJECTORIES 

Anyons are fundamentally defined as two dimensional quantum mechanical systems 
with a specified multi-valuedness for its wave function. The multi-valuedness is stipu- 
lated in terms of one dimensional non-trivial representations of the fundamental group 
of the configuration space, Q. At the least one defines the configuration space for N 
anyons to be, 

Q ~ i?2A. - A, (1) 

where A is the set of points in the plane at which two or more particle positions 
coincide. Denoting by rij the difference of fi and Tj , 

A = {r*j / Tij = for some i ^ j} (2) 

The fundamental group of Q is the so called pure braid group and is known to be 
non-trivial ^. ( For indistinguishable particles Q should be modded by the permu- 
tation group Sn- For most of what follows, this is not essential. ) The multi-valued 
wave functions can be expressed in terms of single valued wave functions by a specific 
explicit multi- valued pre-factor while differentials of multi-valued wave functions can be 
expressed as "covariant" differentials of single valued wave functions. This is discussed 
in detail in reference ^. Using the same notation as ref. 0, in terms of single val- 
ued wave functions, the quantum mechanical Hamiltonian for N-anyons with harmonic 
confinement is given by. 



H = -huiUpl + rl) - «. E ^^ + «^ E "-&} (3) 

^ i i^j 'ij iytj,k 'ij'ik 

— » — * 

Here Pi = — iV. It is convenient to introduce the "statistical" gauge potential Ai 
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A = «,E^x4 (4) 

where k is the unit vector in the z direction. Introducing ttj = pi — Ai , the 
Hamiltonian can be expressed as , 

i 

It is easy to see that the ttj and the r^ satisfy the same commutation relations as Pi 
and fi away from the set A. 

Remarks: 

1) In the above quantum expressions Pi,fi, Uq are all dimensionless. For the classical 
expressions below we will use Uc which has the dimensions of %. These two are related 
by ac = Tiaq. 

2) The use of ttj, -Tj coordinates is the classical counter part of the quantum "anyonic 
gauge" where the Hamiltonian operator expression looks simple and the wave functions 
are multi-valued. 

The corresponding classical system is defined with Q as the configuration space and 
its cotangent bundle as the phase space (F). The Hamiltonian is the same as one given 
in the eq.(3) above. Noting that the Poisson Brackets (PB) among pi^ri are the same 
as those among vfj, -Tj we see that going from pi to ttj is a canonical transformation. The 
classical Hamiltonian may then be taken to be as given in eq.(5). Locally then the anyon 
system is identical to the 2N dimensional isotropic oscillator. The only difference is in 
the global topology of the phase space which turns out to have non-trivial consequence. 

The Lagrangian from which the Hamiltonian, eq.(3), may be obtained is given by, 



^ = ^ E( -1 --1) + «cE% ) ' % = tan-^(^). (6) 

The etc dependent term being a total derivative implies that the equations of motion 
are identical to those of the oscillator. The class of solutions is of course not the same 
due to the removal of the set A. Note also that there is no explicit appearance of ac 
either in the equations of motion or their solutions. Presence of ac only changes the 
classical action and also affects the set of solutions qualitatively. 

The configuration space equations of motion are: 

fi = fi \/i => fij = fij. (7) 

Therefore in the plane defined by any two vectors fi and fj the difference vector 
traces an ellipse in general. For the oscillator degenerate ellipse (straight line) is per- 
mitted but not for the anyons. 

Looking at the Hamiltonian form, one may expect to have the same symmetries for 
anyons as for the oscillator, namely the full 0Sp(4N,R) symmetry. Indeed rank of this 
group being 2N implies existence of 2N constants of motion in involution. These were 
exhibited in reference |^ and were stated to be indicative of the pseudo-integrability 
property. It was pointed out that integrability via action-angle coordinates requires 
further global input. 

III. CLASSIFICATION OF ORBITS AND SYMMETRY REDUCTION 

In this section we will carry out a classification of orbits for A^ anyons and see how 
the pseudo-integrability amounts to a reduction of the dynamical symmetry. We will 
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also see the role played by the fundamental group of the phase space in this regard. 
Noting that a trajectory is specified by giving an initial point, we will use the dynamical 
symmetry of the oscillator system to group various orbits into continuous families. For 
anyons, certain points and hence certain orbits are disallowed (removed) with the result 
that one gets several continuous families. 

A. The case of two anyons 

As usual we can separate the center-of-mass (CM) dynamics from the relative co- 
ordinate dynamics. The CM dynamics is trivial - identical to an oscillator locally and 
globally. The anyonic feature is contained entirely in the relative dynamics described as: 



(8) 



The classical trajectories are trivially known in these cases. These are ellipses with 
two constants of motion, the energy (E) and the angular momentum {J = {r x f)z). 
The sign of J gives the sense of traversal and a trajectory reaches the point r = if 
and only if J = 0. These degenerate trajectories (J = 0) are not allowed for the 
anyons though are allowed for the oscillator. 

Consider the set of all possible trajectories with a given energy. Fix two such trajec- 
tories with angular momenta J, J'. Suppose there exists a continuous family of interpo- 
lating trajectories connecting these two. Clearly if the angular momenta have opposite 
signs, then at least one of the interpolating trajectories must be degenerate. But this 
is disallowed for anyons. Hence, for anyons there must be at least TWO 'orientation' 
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9 = tan^^y/x 


H = 


= i(vr2 + r-) 


7f = p-a^^ 



classes of trajectories distinguished only by the sign of angular momenta. 

While we have shown that there must be at least two families we do not know yet if 
there must be precisely two families. We will use the symmetries to show that there are 
precisely two families for anyons and precisely one family for the oscillator. 

As far as the PB algebra is concerned, the oscillator and anyons are identical and 
hence we do have infinitesimal sjraraetnes forming the Lie algebra of f/(2) (or 05*^(4, R)) 
for both cases. The constant H surface {S^ in the phase space R^ for oscillator) is of 
course invariant under the infinitesimal symmetry transformations. For anyons, the 
phase space has r = removed and the constant H surface has a great circle (one di- 
mensional) removed from the S^, say a "longitude". 

In the appendix A we have given the details of the 0Sp{4N, R) symmetry of the 
oscillator. In the present case of two anyons with center of mass coordinate separated, 
we have effectively, N = 1 and the symmetry generators form 0Sp{4, R). We have 
only the 4 generators of the Tj type. Of these four, ui generates rotations of position 
and momentum while U2 generates time evolution. The generic integral curves are given 
in the appendix A. 

From appendix A we also recall that 0Sp{4, R) acts transitively on H = E sphere. 
This shows that for the oscillator there is precisely ONE family of trajectories. For 
anyons we need analogous result with the extra condition that degenerate trajectories 
are avoided to conclude that there are precisely TWO families. 
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Claim Given any two trajectories with the same sign of their angular momenta, there 
exist a one parameter group connecting the two without changing the sign of angular 
momenta. 

Proof Observe that ui,U2 transformations leave the angular momentum invariant. 
Choose u = cos(/3)u3 + sin{j3)u4. The one parameter group generated by this 
generator transforms the angular momentum as, (appendix A) 

2J{a) = 2J - sin^{a) {2J - ujuLvluj} + sin{2a) {cuLucu} (9) 

For the choice of u made, uLu = — L and putting J„ = "^^^^ ^ 

J{a) = J cos{2cr) + sin{2(r) J„ 

= \/j^ + Ju cos{2a - 6) , S = tan'\Ju/J) G (-7r/2,7r/2) (10) 

This is valid for all choices of 13. Observe that J„ depends on (3. Furthermore it 

satisfies, 

d^Ju{l3) 



dp^ 



- Jui(3) 



(97 

Ju{f3) = Ju{0)cos{f3) + ^{0)stnif3) => (11) 

The last equality follows by explicit evaluation. Choosing (3 to maximise J„ then 
implies that the square root in the equation for J{cr) is just the energy E. Since for 
any given energy, we must have J^ < E'^, it follows that J(o') covers the entire range of 
possible J values. 
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From this it follows that we can always choose a a such that J((t) = J' without 
passing through zero angular momentum. This proves the claim. Hence there are pre- 
cisely TWO families of trajectories each spanning the whole orientation class. This gives 
the classification of trajectories. 

One may view the S^ as compactification of R^ with the south pole as the origin of 
R^ and the north pole as the point at infinity. The removal of A now corresponds to 
removal of an infinite line, say, the z-axis of the R^. It is immediate that the fundamental 
group is the group of integers, Z, and therefore the basic trajectories winding around 
the z-axis belong to two orientation classes. In this case the fundamental group directly 
shows how the basic trajectories are naturally classified into the "orientation classes" . 

Consider now the Hamiltonian vector fields generating the three one parameter sub- 
groups of S'f/(2),ui generates rotations about the z-axis. Its orbit through any generic 
point avoids the z-axis and the vector field remains complete (its integral curves range 
over the full real line) . The other two subgroups on the other hand have orbits necessar- 
ily cutting through the z-axis. These vector fields therefore are necessarily incomplete 
(for the anyon case). The transformations generated by these vector fields do not 
for a group and the dynamical symmetry group for anyons is reduced from U{2) to 
f/(l) X [/(I). 

Thus the removal of the set A has altered the fundamental group, has classified tra- 
jectories into orientation classes and has also reduced the dynamical symmetry group 
form f/(2) to U{1) x f/(l) due to the incompleteness of vector fields. 
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Note that while the symmetry group is reduces from f/(2) to f/(l) x [/(I) for anyons, 
the rank has remained the same and therefore the anyonic system continues to be inte- 
grable. The integral curves being circles, we have integrability via action-angle variables 
and we get the exact spectrum from the semiclassical approximation. In fact one can 



follows the EBK quantization procedure to reproduce the exact spectrum. ||T3| 



B. The case of many anyons 

The N > 2 case differs from the N = 2 case in many ways. The potential symme- 
try group, 0SP{4:N, R), (see appendix A) generating families of trajectories is lot more 
complicated and so is the fundamental group for the phase space of anyons. However the 
analysis of the previous section already suggests a strategy for obtaining a classification 
of the families of trajectories. A generic trajectory may be viewed as a collection of 
^ ^~ ' ellipses traced by the fij,i < j. To get a convenient handle on the disallowed 
trajectories, define Jij = {fij x ry)^. Clearly a trajectory will cut through the set A if 
and only if J^j vanishes for at least one pair of indices. We will term such a trajectory 
as degenerate. For anyon only non-degenerate trajectories are allowed. Note that Jij 
are constants of motion just as Jj and Ei = ^{rf + p'f) are. Each of these elliptical 
trajectories will have a sense of traversal ( e^j = sign of Jij ). Thus a basic trajectory of 
A^ particles has an associated set: {e^}. Exactly as in the previous sub-section, the set 
of non-degenerate trajectories can be classified by these 'orientations'. Provided we can 
exhibit trajectories which realise all possible choices of eij, we will have at least 2^(^~^)/2 
families for anyons and precisely as many if every member of a class is connected to ev- 
ery other member of the same class by a symmetry transformation without crossing A. 
It is easy to show that every choice of eij is realised by considering concentric circular 
trajectories for each of the particles with various possible signs for the angular momenta 
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Ji and various possible ordering of the radii of the circular trajectories. For oscillator 
we will have precisely one family if 0Sp{4N, R) acts transitively on the constant energy 
sphere, S*'^~^. 

The proof of transitivity in the case of oscillator is most directly obtained by noting 
that 0Sp{4N, R) is isomorphic to U{2N, C) and it is well known (and easy to see) that 
the coset space U{2N,C)/U{2N — 1,C) is diffeomorphic to S"^^^^. This immediately 
shows that for the oscillator, there is precisely ONE family of periodic trajectories. The 
counting of families is more tricky in the case of the anyons. Let us see the symmetry 
reduction in a slightly different manner. 

Let A' denote the set of all degenerate trajectories i.e. at least one of the Jjj's is 
zero. Note that A and A' are different, the former being a proper subset of the latter. 
Suppose we find the subgroup of symmetries which leave A' invariant, then the same 
subgroup will also leave the set of all non-degenerate trajectories invariant. Since we 
are looking for a group, we can consider the infinitesimal action. Let J = Hi^j Jij- 
Clearly, J = characterises the set A'. If we have a degenerate trajectory with two 
or more J^'s zero, then every infinitesimal action will keep within the set of degenerate 
trajectories. If only one Jmn is zero then infinitesimal action will act invariantly if and 
only if Jmn is invariant. Of course we could have degenerate trajectories with different 
angular momenta being zero and hence for the invariance of A' it is necessary that the 
generators of the subgroup must leave each of the Jij 's invariant whenever Jij itself is 
zero. This subgroup is identified in the appendix A. There are two forms of generators. 
Any generator with all the block matrices being the same u matrix with u G 05*^(4, R) 
will leave all the J's invariant. These transformations act on the center-of-mass variable 
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alone and constitute the expected 0Sp{4, R) symmetry. In addition, generators with 
only diagonal blocks being the same u and all the off-diagonal blocks being also leaves 
the J's invariant provided u = Ui or U2 (or a linear combination thereof). These are 
just the total Hamiltonian and the total angular momentum {J2i Ji) which generate time 
evolutions and common rotations of all the positions and momenta. These are the only 
subgroups leaving the set of degenerate (and non-degenerate) trajectories invariant. The 
vector fields corresponding only to these will remain complete. Thus we see that the 
symmetry is reduced to 05*^(4, R) x 0(2) x 0(2). This is identical to the symmetry in 
the case of two anyons. 

It is apparant from the above discussion that pseudo-integrability of anyons really 
means that not all the Lie algebra level symmetries exponetiate to Lie group symmetries. 
One should distinguish between "infintesimal integrability" and integrability. A general 
discussion is given in the last section. 

This purely classical analysis has already explained the qualitative result that for 
anyons only the total energy and the total angular momentum are the conserved quanti- 
ties (good quantum numbers). Since we also have a classification of periodic trajectories 
we will take the next step of attempting semiclassical approximation to see what further 
information can be obtained. 



IV. SEMICLASSICAL SPECTRUM FOR ANYONS 



The central quantity of interest is the propagator defined as, 
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G{E + te)^Tr{ ^ 



E-H + ie' 

= Y.{E-En + ie)-^ (12) 

n 

= J2niE- E„)-i) - m 5] 6iE - E„) 

n n 

In the last equality, V denotes the principle value while the second term is of course 
the density of states. Defining G = iG allows us to write the density of states, g{E) as , 

g{E) = ^{G{E + ze) + G*{E-ze)} 

= - Re{ G{E + le) ) (13) 

TT 

The derivation of the semiclassical trace formula begins with the definition of the 
propagator, G{E), with trace operation expressed via a path integral representation. 
The trace operation combined with the stationary phase approximation (SPA) gives the 
propagator as a sum over periodic orbits (in the phase space), of terms with an amplitude 
given by a suitable Van-Vleck determinant times a phase whose exponent is the value 
of the 'action', / p • dq around the periodic orbit together with a contribution from the 
Maslov indices. This works very well for isolated periodic orbits. However for systems 
with symmetries, as for the oscillator or anyons, the periodic orbits come in continuous 
families and a modification is needed. 

In such cases, the sum over orbits gets replaced by a sum over families of orbits 
together with a measure factor. If the family arises due to a symmetry, as is usually 
the case, the amplitude and the phase is same for all members of the family and as 



such can be computed from any one member. This is discussed by Littlejohn et al |10 
in detail. The origin of families of periodic orbits in the case of oscillator is of course 
the 0Sp{4:, R) dynamical symmetry and for anyons essentially the same transformations 
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generate families of trajectories. This is already discussed in the previous section. 

To obtain the action along a periodic orbit one point needs to be noted. Observe 
that in terms of the 7fj,ri coordinates there is no explicit reference to a and A is also 
precisely the set on which one gets (5-function singular PBs with a appearing as the 
coefficient. If one did inverse Legendre transformation to get a Lagrangian one would 
not get the a dependent term, but of course the configuration space will have the A 
removed. Using pj, -Tj variables and then doing inverse Legendre transformation will pro- 
duce the Lagrangian with the a dependent total derivative term. Now the removal of 
A is made explicit by the a dependent term. It is this explicit form which is convenient 
for computation of action for periodic orbits. 

The action for the oscillator is equal to {2txE)/uj and the action for anyons is the 
same as that for the oscillator except for the contribution from the total derivative term 
in the Lagrangian. Its value is ±27ra with sign depending on the orientation class(es). 
We have already seen that for A^ = 2 we have precisely 2 orientation classes while for 
general A^ we have at least 2^(^~^)/2 classes. 

Observe that even without the knowledge of precise number of families, we have 
precisely 2^(^~^)/^ sums in the sum over periodic trajectories. This is because, the con- 
tribution of the action integral depends only on the orientation class and not on further 
possible subclasses of trajectories. The measure factors therefore add up within each 
orientation class and this is sufficient to get the semiclassical eigenvalues. 

The sum over families in the trace formula contains as many terms. To get the eigen- 
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values we need to look at only the action integral which differs from the oscillator only 
by the additional / ^ ^i^j ^ij contribution. The oscillator's contribution to the action 
integral is of course just —E, where E is the classical energy and u is the oscillator 
frequency. 

The net ac dependent contribution to the action integral is given by, 

|^/e% = 2™, 5: 6., (14) 

Here eij is the sign of traversal of r^j and is ±1. The J2i<j is trivial to evaluate. 

Since the classical trajectories are identical to those of the oscillator, we will have the 
same Maslov indices which incorporate the usual zero point energy namely, A^(A^ — 1) 
for 2N dimensional oscillator. Including this contribution, the semiclassical energies are 
given by : 

En{ei} = hu{n + a,Y.e^3 + N{N-l) ) n >0 (15) 

The coefficient of Uq takes all possible integer values from —N{N — l)/2 to 
A^(A^ — l)/2 in steps of 2. The extreme values correspond to all the particles traversing 
the same way and represent the contribution of "collective" motion shown in 0] . The re- 
maining values represent contributions of "relative" motion. These are new eigenvalues 
albeit at the leading h level. Numerical spectra available for A^ = 3, 4 show eigenvalues 
matching with the above at the a, = and a, = 1 0]- These numerical eigenvalues 
however are non-linearly interpolating while our semiclassical eigenvalues are linearly 
interpolating. It is possible that the semiclassical eigenvalues will receive corrections 
to make them non-linearly interpolating. In reference P| also it was conjectured that 
there will be linearly interpolating eigenvalues with various slopes based on a different 
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semi classical argument. We have obtained these eigenvalues by direct application of the 
periodic orbit theory and proved their existence for all N. 

Observe that the aq dependence in the eigenvalues ( locations of poles of G{E) ) 
arises only from the explicit ac = cxqh in the action. This therefore must be at the most 
linear in Uq. Since Uq is fixed, the semiclassical limit of ^ — > implies etc ~* 0. Thus 
one may at the most see linear Uc dependence in the leading approximation which could 
however be indicative of exact non-linearly interpolating eigenvalues. 

Thus the leading level semiclassical propagator already indicates the presence of 
eigenvalues such that E{aq = 1) — E{aq = 0) takes all possible integer values between 
±A^(A^ - l)/2 in steps of 2. 



V. HALF PERIOD TRAJECTORIES AND EXACT PROPAGATOR FOR 

TWO ANYONS. 

Recall that anyons are fundamentally defined quantum mechanically and all the en- 
ergy eigenfunctions of the system vanish on the set A of coincident points. We modeled 
the classical system by removing A and explored the consequences in the previous sec- 
tions. In particular we just omitted trajectories that could cut through A. 

From a purely classical point of view, though, this is little unsatisfactory. One could 
legitimately ask just what happens to trajectories (eg J = ioi N = 2) that attempt 
reaching the disallowed region? To decide this one has to extend the classical modeling 
by supplementing the equations of motions with a specified "boundary" condition. The 
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choice is helped by noting that the curl of the vector potential is a (5-function which may 
be "regulated" to stipulate the "boundary condition" . This analysis is given in appendix 
B. 

The result is that all the elliptical trajectories are identical to those of the oscillator 
and only the degenerate trajectories get modified. A degenerate trajectory must reflect 
back from the coincident point. 

Clearly such an trajectory will have half the period of the generic trajectories and 
hence we will refer to these as the "half trajectories" . 

Should these orbits be included in the sum over orbits in the context of semiclassical 
approximation? Are they already "visible" in the exact spectrum known for two anyons? 
To answer these, let us look at the exact propagator for the (relative coordinate part of) 
the two anyon system. 



The propagator can be computed exactly from the known exact spectrum [jrT| . The 
exact spectrum is given by, 

En,j = huj{2n + \j - aq\ + 1), n > 0, j eZ (16) 

From this the partition function is obtained as, 

v(R\ _ '^osh{ f3huj {aq - 1) ) + cosh{ (3hujaq ) 

^^^^ - 2stnh\phco) ^^^> 

and the exact density of states is obtained as: 
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E , ,^, .2TxkE , , .2TxkE , ,^ 

--— -^{1 + > ,|cos(— ^ h 2Txkaq) + cos[-— 2Txkaq)\ 

(nujY ^ nuj nuj 

1 2'KkE 
— —_ — - y^|2ao sin(2Trkaa) sin(^ )| (18) 

^{(—1)^ sin{Tikaq) sin{— — )} 



(huj) ^^ huj 

Rewriting the products of sines as differences of cosines we get, 

E 



9{E) 



I , E ,^ ,2'KkE , , 

+ —- ( ^ + ttg ) V cos[-— h 2Txkaq) 

nuj nuo ^ nuo 

I , E ,^ ,2'KkE , , 

+ —- ( ttg ) V cos[^ 2'nkaq) (19) 

Remarks: 

1. The partition function is manifestly invariant under Uq -^ l—Uq . The convergence 
of integrals in computing the density of states for Uq G (0, 1) requires ^ to be greater 
than or equal to 1. 



2. The first term above is the usual Thomas- Fermi term ||TT|. This will be suppressed 
in expressions below. 

3. In the limit /i — > keeping E, Uq fixed, the exact propagator effectively looses all 
dependence on aq. This is also the semi classical trace formula for the oscillator showing 
that the semiclassical approximation is exact for the oscillator. All the other terms are 
sub-leading. 

4. In the semiclassical computation, one will use Uc = haq and the SPA will be done 
keeping etc fixed. Then only the last two terms will be sub-leading and the propagator 
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will continue to have Oc dependence. 

5. In the last two sub-leading terms the argument of cosine is half of that in second 
and the third term. As is well known, the argument of cosines is the classical action, 
J p- df, around a classical orbit including possible contributions from "Maslov indices". 
The last two terms are therefore suggestive of contribution from "half orbits" . Indeed 
if one takes a Fourier transform of g{E) []TT| , then one sees a peak at 1/2 period from 
the last two terms. One may therefore conclude that half orbits should also be included 
in the trace formula. 

This turns out to be somewhat ambiguous. To see this let us express the propagator 
in a sum-over-orbits form. To do this rewrite the the last two terms to resemble the 
second and the third terms. 

The last two terms also have a (—1)'^ in the summation over k and this can be han- 
dled in two ways. 

(A) Separate these sums into even and odd integer sums. All the sums over k being 
geometric series can be done explicitly (ie needs to be added). Defining S = E/{huj), 
we get, 

hujg{E) = S + 



1 — Q2mag{£+aq) 

- (£ — a + - \e'^M£-aq) _ l^in{£-aq) 

]_ Q2niaq{£~aq) 

Here C. C means complex conjugation. Now we can "read-off' G{E + ie) and get 
(suppressing the Thomas- Fermi term), 
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1 1 p2TTi(£+aq+ie) 

1 1 p2TTi{£—aq+ie) 

2 i^ - «^ + 2^1 - e-(^-^+-)) } { ^_^,..(,_,^^,) } (21) 
Noting that the second group of braces is sum of a geometric series: 

„2ni{£±aq+ie) 



1 p27ri(f ztog+ie) 

-^ ^ k>l 

we see that the propagator is expressed in a form suggestive of a sum over periodic 
orbits. The first group within the braces being the "amphtude" while the second group 
being contributions from muhiple traversals of a basic periodic orbit. The two terms can 
be seen to come from famihes of eUiptical orbits going clockwise and anti-clockwise sense 
and that there no term corresponding to degenerate ellipses or half orbits. Note that 
the terms in the "amplitudes" other than S are sub- leading relative to S. This leading 
part of the "amplitudes" can also be seen to come from POT in presence of continuous 
families of periodic orbits [|T0|. The "amplitudes" however are complex. For a^ = 
there is cancellation of these sub-leading terms and one recovers the exact result for the 
oscillator. It appears that though the system is integrable, the trace formula does not 
give exact propagator. 

The poles in G{E) come from the two terms and are given by: 

i?+(r7,_|_) = huj{n+ — Oq) ^+ > 2 (23) 

E_{n_) = huj{n^ + aq) n_ > 1 (24) 

The residues at these poles { £ ±aq = n± ) are n± if n^ is even and n-t — 1 if n^ is 
odd. These residues of course give the degeneracy. Note that the sub-leading terms are 
important for this. 
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It is easy to see that the locations of all the poles can be re-expressed in the form 
given by the exact spectrum. The degeneracies at the poles also match exactly as ex- 
pected. 

(B) We can also write {—l)''cos{k9) = cos{k{6 ± tt)). Proceeding exactly as in the 
case (A), we can "read-off" G{E + ie) and get, 



-huj\ 



2711 I 



1 „2TTi(£+aq+ie) i p2m{£—aq+ie) 



1 ^ni{£-l+aq+ie) 



-7{ 



4 '- ]_ QTri{£-l+aq+ie) 



} + tI 



1 pTri{£+l-aq+ie) 



4 ^ X — QTri{£+l-aq+ie) 



} 



(25) 



These have poles which are subsumed by the poles from the first two terms. The 
overall pole structure and residues of course are exactly same as before as they should 
be. 

The exponents in the last two terms however do not look like the action integral. A 
reflecting trajectory will not receive a contribution from the ag dependent total deriva- 
tive term. The reflecting trajectories will also form a single family of trajectories and 
thus should give only one term and not two terms. The last two terms can be combined 
and the propagator may be re-expressed as, 



-^—\ G(E + ie) = y i ^ + '^q \ ^27Tik{£+aq+ie) _^ /^ ^x ^2nik{£-aq+ie) 



-isin{kTxa) i^k{£~i) 
^ 2 ^ 



(26) 
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The propagator now does have form indicating a contribution from the half period 
orbits including a contribution from a Maslov index due to reflection. However it has a 
complex amplitude with a dependence on a^ and the multiple traversals index k. 

Thus the exact propagator can be expressed in two different forms mimicking the 
trace formula with and without the half period trajectories. Although the phases are as 
expected, the amplitudes are not. These terms also have different orders of Ti. 



The different orders of fi seen can be understood in the context of trace formula as 
due to families of periodic trajectories and not due to any higher order corrections to the 
SPA. When periodic trajectories come in continuous families, the number of Gaussian 
integrations is reduced since the integration over the family must be done exactly [10 



and this increases the powers of Ti in the denominator. The elliptic trajectories form a 
three parameter family while the half trajectories form a one parameter family which 
explains the difference in the powers of Ti. 

While the powers of h, can be understood, the precise matching of the amplitudes 
from the POT does not seem to follow. The appearance of complex amplitudes may be 
indicative of the need for complex trajectories and/or precise matching of semiclassical 
wave functions across the excluded regions. This is beyond the scope of present work. 

The ambiguous "presence" of half trajectories seems to suggest that one may either 
exclude these altogether but extend the trace formula to deal with the non-trivial topol- 
ogy due to the excluded regions or include these trajectories as additional families but 
perhaps use complex trajectories. Needless to say that the structure of half trajectories 
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will be more complicated for N > 2. For getting the semiclassical eigenvalues, the first 
alternative may suffice and is already adequate to indicate the presence of non-exactly 
known eigenvalues. As there does not appear to be any scope for going beyond the 
leading SPA (all intermediate exponents of phases are quadratic), it is not clear how one 
may develop a semiclassical scheme for obtaining the non-linear a^ dependence. 



VI. SUMMARY AND CONCLUSIONS 

In the previous sections we encountered two important features, one related to the 
fundamental group and one related to symmetry reduction due to incomplete vector 
field. Both were caused by the same source, namely removal of coincident points imply- 
ing topologically non-trivial phase space. Both are relevant for periodic orbit theory in 
a general way and a few general remarks are in order. 

The stationary phase approximation to the propagator naturally leads to periodic 
orbits in the phase space. These orbits may be isolated or come in continuous families 
or both. The continuous families may be generated by a full group of symmetries or by 
only a subset of symmetry transformations. Since each periodic orbit is also a map of 
S^ to the phase space F, clearly every orbit must belong to one and only one homotopy 
class of the fundamental group, vri(F), of the phase space. If an orbit is a member of a 
continuous family then the entire family must belong to a single homotopy class. Note 
that a given homotopy class may contain no orbit (solution of equation of motion) or 
several isolated orbits and/or several families of orbits. But a family can not spill over 
two distinct homotopy classes. In section 111 we saw precisely the splitting of a single 
basic family for oscillator into two basic families for anyons because of the non-trivial 
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fundamental group. Multiple traversals of course belong to different homotopy classes 
and are explicitly summed over. Thus a non-trivial tti may (but not necessarily) provide 
an obstruction to a symmetry. How exactly may such an obstruction manifest itself? 
For this we have to consider vector fields generating symmetries. 

Recall [jT^ that every function on the phase space generates infinitesimal symplectic 
diffeomorphisms (canonical transformations) via its corresponding (globally) Hamilto- 
nian vector field. The Lie algebra of such vector field is isomorphic to the Possion 
Bracket (PB) algebra of functions on F. However such infinitesimal transformations 
exponentiate to give a one parameter group of transformations only if the vector field is 
complete i.e. the integral curves of the vector field can be extended so as to have these 
as a map from the full R. Only complete vector fields - and this is a global statement 
- give rise to groups of symmetries. (A corresponding quantum mechanical statement 
for continuous symmetries is contained in the Stone's theorem: every one parameter 
group of unitary transformations is generated by a self-adjoint operator and conversely.) 
Since we usually want classical symmetries to be reflected at the quantum level with 
observable generators we have to have groups of symmetries and hence complete vector 
fields. 

The criterion of integrability in terms of vanishing PB's guarantees only the exis- 
tence of infinitesim,al sym,m,etries which is of course a prerequisite. An "infinitesimally 
integrable" system may thus be: (a) integrable via action-angle coordinates if the vec- 
tor fields are complete and the integral curves are periodic; (b) integrable, but not by 
action-angle variables if the vector fields are complete but only a subset of these have 
periodic integral curves and (c) partially integrable if only a subset of vector fields are 
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complete. The "pseudo integrability" property of anyons pointed out in ref. [Q falls in 
the category (c). As explained in the appendix A apart from the Hamiltonian only the 
total angular momentum has a complete vector field and hence is the only quantum 
number that survives. 

A non-trivial fundamental group by itself however does not imply possibility of in- 
complete vector fields. For, on a compact manifold all vector fields are complete and it 
can of course have a non-trivial tti. In many cases constant energy surfaces are compact 
and one does not have to worry about incomplete vector fields. For the anyons however 
removal of the set A makes the constant energy surface noncompact which admits pos- 
sibility of incomplete vector fields and corresponding loss of symmetry. 

In the present work we have seen manifestation of all these features. 

We also considered application of the periodic orbit theory to this locally trivial 
(integrable) but globally non-trivial system. As a by-product, we saw how the exact 
spectrum for two anyons is reproduced by semiclassical methods. We saw that while the 
dynamical symmetry is reduced from SU{2) to f/(l), the rank remained the same and 
hence integrability property is preserved. 

For A^ > 3, we reproduced the previously known exact eigenvalues [^]. In addition, 
we obtained further new linearly interpolating eigenvalues. In the language of reference 
Q, these correspond to effect of 'relative' dynamics. In this case the symmetry reduction 
was drastic, from OSp{AN, R) to OSp{A, R) x 0(2, R) x 0(2, R). The rank was reduced 
from 2A^ to 6, destroying the integrability property. 
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We discussed in detail the issue of "half trajectories and pointed out their am- 
biguous role. Taking the view-point that "half trajectories be excluded, and noting 
that there does not seem to be any scope for computing higher order corrections to the 
semi classical spectrum with a possible non-linear dependence on aq, it seems that the 
non-linearly interpolating eigenvalues are genuine quantum consequences beyond what 
semiclassical analysis could give. Semi-classical analysis is nevertheless sufficient to in- 
dicate the presence of these eigenvalues. 
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Appendix A : OSp dynamical symmetry and classification of trajectories 

In this appendix we collect together some of the well known facts about the dynam- 
ical symmetry of the n dimensional isotropic oscillator. 

The isotropic oscillator in n dimensions has 0Sp{2n, R) as the group of dynamical 
symmetries. This is a group of 2n x 2n order matrices, g, which are both orthogonal 
and symplectic. Denoting by g the transpose of g, we have the defining equations: 

g g = hn ■ (orthogonal) 

g 17 g = 17 : (symplectic) 

and with g ^ l2n + eT generators T satisfy (27) 

T = -T 

T 17 = -17 T 

Here l2n is the identity matrix of order 2n while 17 is a suitable matrix defining the 
symplectic condition. This will be chosen below. 

The symplectic condition ensures that we have a (linear) canonical transformation 
while orthogonality ensures that the Hamiltonian, J^ApI + (li)/'^ , is invariant. It is easy 
to see that dimension of this group is n'^ while its rank is n. In fact one can show that 
the group 0Sp{2n, K) is isomorphic to the group U{n, C). Making an explicit choice of 
fi in a block form, one can obtain the block form for g using the symplectic condition. 
The real block matrices can be combined into complex block matrices. The orthogo- 
nality condition in terms of real matrices then translates into the unitarity condition 
for the complex matrices thereby proving the group isomorphism. This isomorphism 
immediately implies that the ortho-symplectic group acts transitively on the constant 
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energy sphere. This is used in the section III. 

We are interested in finding the action of this group on the phase space via canonical 
transformations (symplectic diffeomorphisms) . For notational convenience let us group 
the standard canonical variables together and denote them by o;^ 1 1^ = 1,2, ...,2n, the 
first n being coordinates and the last n being the momenta. The dual of the sym- 
plectic form has components, Vt^^ given by the ^v th element of the matrix f2. With 
this notation, the PB of functions on the phase space and the infinitesimal canonical 
transformations are given by, 

{F{u) , G{uj)} = Q'^' d,F d,G 

6,uj^ = e n^"" d^F{uj) (28) 

Functions purely quadratic in u generate linear canonical transformations and are 

also closed under the PB's and those which leave the Hamiltonian invariant give the 

symmetry transformations. Explicitly, 

F = ^cu A u , G = ^cu B cu =^ 

' ' (29) 

{F, G} = \u (AflB - BflA) u 
The Hamiltonian corresponds to A = l2n- The Poisson bracket of any G with H van- 
ishes provided the matrix B commutes with f2. Thus generic matrices defining quadratic 
functions are real, symmetric matrices commuting with f2. The matrices f2B are then 
antisymmetric and provide an isomorphism of quadratic functions to the generators T 
of the group 0Sp{2n, R). 

Action of the one parameter group generated by a function G is found from the 
integral curves of the corresponding Hamiltonian vector field. These curves are defined 
by the matrix equations, (T = ilB) 
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^ = flBc (30) 

da 

uj{a) = (e'^"^)cc)(0) (31) 

It is convenient to choose a particular grouping of the phase space coordinates and 
corresponding choice of fi. Firstly let us put n = 2N, relevant for the present context, 
so that the phase space is 4N dimensional. Arrange the coordinates and momenta as 
Xi,yi,pij:,piy, ....,Xiy,yN,PNx,PNy and denote by cJj the coordinates and momenta of 
the i^^ particle. The index i now runs from 1, ..., N and each o^j is a 4 x 1 matrix. 
Correspondingly we choose f2 as a block diagonal matrix with N blocks and each block 
being 4x4 matrix A. We choose, 



/ 



A 



O2 l2^ 



(32) 



y -I2 O2 J 



One can choose a basis for 0Sp{4N, K) as follows. Let us denote the generators as 
Tj, i = 1,2, ..., A^ and Tj^ with i < j. Each of these are expressed in the block form. 
The Tj are block diagonal with a nonzero 4x4 matrix Uj as the i*^ block element. The 
Tij have a matrix Vjj at the i*'' row and j*^ column {{ijY^ block) and — Vjj at the {jiY^ 
block. In equations. 



\-^i)mn Uj OifYi Omn 

With these definitions it is easy to translate the conditions on the generators in terms 
of the 4 matrices u , v as: 

Ui = -Ui UijA = Aujj ; Vjj A = Av^j (34) 

Thus the u's generate 0Sp{4, R) while the v's are required to commute with A. The 
number of independent u's is 4 while number of independent v's is 8 of which 4 are 
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'diagonal' and 4 are 'off-diagonal'. The dimension of 0Sp{4:N, R) is thus 4A^ + 8A^(A^ — 
l)/2 = 4N'^. These independent matrices can be explicitly chosen in terms of 2 x 2 
Pauli matrices and the identity matrix as: 

/ ia2 O2 1 / O2 I2 1 (02^1) [ O2 ^3 1 



U(l,2,3,4) 



V(l,2,3,4) 



V(5,6,7,8 



ia2 U2 
O2 icr2 



V 



-'2 -12 
I2 O2 



V 



h 0"! 
-O"! O2 



— 0-3 O2 



/ 



/^ I2 O2 ^ ( a^ oA I ido 02^ ( (73 O2 ^ 



12 W2 
O2 I2 



O"! U2 
O2 0-1 



O2 io-2 



/^ O2 I2 ^ /^ O2 (71 ^ I Q 



>2 A2 
-I2 O2 



2 O"! 
(71 O2 



-I2 «(7'2 
-i(72 O2 



y O2 (73 ^ 
^ ^02 .3^ 



(35) 



V 



h (73 

-(73 O2 



(36) 



We are interested in one parameter groups generated by some element T of the Lie 
algebra. If T satisfies: T^ = -P with P satisfying P^ = P , [T, P] = 0, then it follows 
that, 



^o-T 



I4JV - P + P [cos ((7) + sm((7)T] P 



(37) 



For the basis generators Tj we have P to be block diagonal with I4 as the i^^ block 
while for Tjj we have P to be the block diagonal matrix with I4 as the i^^ and the f^ 
blocks. Using these, the exponentials of basis generators can be evaluated to get the 
integral curves as, (with obvious notation) 

oj{a) = [!'■ + li { cos{a) + sin{a) T^ } Ij ]cj(0) and, 

uj{a) = [ I'ij + lij { cos{a) + sin{(T) Tij } lij ] u{0) (38) 

Therefore integral curves of the basis generators are periodic curves. Further, the 
Tj's affect only the i^^ particle position and momenta while the Tjj mix the i*'' and j^^ 
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particles only. 

We need to study how the angular momenta Jj and Jij vary under the action of one 
parameter subgroups. Define the (block) matrices, 

\^ij )mn = ^ \ ^mny^irn ~r Ojm) OiftiOjn OinOjm f V'^'^J 

In terms of these matrices the angular momenta are given by: 

1 _ 1 _ 

Ji = - u Li uj = - uJiL Ui 

1 _ 1 _ 

Jij = - u Lij u = - LOij L Uij , Uij = Ui -Uj. (40) 

Under the group generated by the basis generator, Tj, the Jm for instance varies as: 

2Jm((T) = 2Jm + 5im —sin'^{(j) {2Jj — a); (u L u) uj{\ 

+ sin{a)cos{a) { u-^ (uL + Lu) cj; } ] (41) 

while under the group generated by the basis generator, Tj^ the Jm varies as: 



'^•JfriK,^} '^•Jm ~r Oi^yi 



(jjrn 



—sin (a) {2Jj — cjj (v L v) uj-j] 

+ 2sin{a)cos{a) { uJi (Lv) ujj } ] 

—sin^{a) {2Jj — oJi (v L v) cui} 

- 2sin{(j)cos{a) { Ui (vL) cjj } ] (42) 

Similar but more complicated expressions follow for the Jmn{^) also. 

Remark: For the case of A^ = 1 we have only one block. This could be either a 
single two dimensional oscillator OR the relative coordinate dynamics of two anyons. 
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The generators are of course only u's. Of these U2 corresponds to the Hamiltonian itself 
while ui generates same rotations of both r*, p. These two matrices commute with L 
while the remaining two anti-commute with L. This leads to the result quoted in the 
section III . 

To deduce the surviving symmetry group for many anyons we need infinitesimal 
variations of Jmn for arbitrary generator T. This is easily derived and is given by, 

^Jmn = 2^ ^mn L (T^^ — T„j )lOj (43) 

3 

In section III we needed to determine T such that 5 Jmn is zero whenever Jmn is 
zero. That is, at all points in the phase space where any single Jmn is zero, we want its 
infinitesimal variation induced by T to be zero. 

Fix a particular m and n. Fix oOmn- We can consider points with u;^, j ^ m,n such 
that no other Jij is zero. SJmn = then implies Tmj = Tnj V j 7^ m, n. We can 
also consider points with different Um, uJn keeping Umn fixed and maintaining all other 
conditions. This implies Tmm ~ T„m + Tmn ~ T„„ = 0. Repeating this for all m ^ n 
fixes the form of T in terms of arbitrary generators of 05*^(4, i?), u , v as, 

{T)ij = u5ij + v(l-%) , [L, u-v] = (44) 

The choice u = v corresponds to Tjj = u V i,j. It follows that all Jmn 's are 
invariant independent of their values. It is easy to see that these four T's effect trans- 
formations of the center-of-mass variables which are insensitive to the anyonic features. 
This 0Sp{4, R) symmetry is thus always present for all N > 2. 
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The choice v = imphes T is block diagonal with the same u on all the diagonal 
blocks. Further, [L, u] = implies that u must be a linear combination of ui and U2 . 
These two T 's can be seen to correspond to the total Hamiltonian and the total angular 
momentum. This result is used in the section III to deduce that the surviving symmetry 
for N anyons is 0Sp{4:, R) x 0(2, R) x 0(2, R). 
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Appendix B : Regularised classical dynamics and reflecting orbits 

Consider without loss of generality the case of two anyon in the relative coordinates. 
Generic orbits in the configure space, R^ — {0}, are of course an ellipse and only a de- 
generate elliptical orbit attempts to pass through the origin. Noting that this system 
can be thought of as a charged particle in presence of a singular ( "statistical" ) magnetic 
field at the origin, one may regularise the magnetic field or the flux to study the orbits 
and obtain the limiting behaviour to deduce "boundary condition" at the origin. 

To do this we imagine the relative Hamiltonian to be that of a particle in a magnetic 
field along the z axis. The field is of course to be effectively confined to a small disc 
around the origin. Now observe that for an axially symmetric magnetic field along the 
z-axis, B{r, 9) = B{r)k where r, 9 are the usual spherical coordinates in two dimensions, 
we may write the vector potential as, 

A = Arf + Ae9 (45) 



This implies. 



B{r) = drAe + ^ - ^9^^,, (46) 



In a symmetric gauge, the vector potential is independent of the polar angle and by 
choosing it to be divergence free one can set the radial component to zero. 

The flux $(-R), through a disc of radius R is given by, 

<^{R) = 2tt f drrB{r) = R i Aed9 = 2iiRAe{R). (47) 

This implies, 

37 



Ae{r) = ^ = lJ^dr'r'B{r') (48) 

The choice, $(r) = 27ra V r > 0, gives the vector potential used in the quantum 
mechanical calculation. It also implies the rB{r) = a6{r) and this of course is the 
singular nature of the magnetic field. Notice that smearing the 6{r) will not make the 
magnetic field non-singular at the origin because of the explicit 1/r. 

For a regulated system one wants the magnetic field to be non-singular every where 
and effectively confined to a disk of radius e around the origin. This requires the 
vector potential Ag also to be non-singular and therefore vanishing at the origin. If 
Ag{r) -^ cr^ as r approaches zero then, B{r) -^ c{(3 + l)r^^^. For a finite, nonzero 
B{0) we must have (3 = 1 and therefore <I>(r) -^ 27rcr^ near the origin. For r > e we 
still retain the flux to be 27ra. Continuity at r = e then gives c = a/e^. Notice that this 
limiting behaviour is fixed by the demand of non-singularity of the fields and as such 
must be reflected in any explicit choice for the magnetic field. 

Thus our regulation involves choosing. 



$(r) = 27ra V r > e 

27ra^ as r ^0 (49) 



e^ 



One could choose a uniform nonzero magnetic field inside the disk as an explicit 
choice but it will not be necessary. Since we are interested in the limiting behaviour of 
trajectories as e is taken to zero, the limiting behaviour of the flux is all that we need. 

Consider now the orbit equation for e nonzero. Integrating the equations of motion 
once using the two constants, energy E and angular momentum i, the orbit equations 
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in r, 9 coordinates become : 



r 


= 


± 


^2E-7 


-.2Q2 


_^2 


e 


= 


i- 


- rAe{r) 

,2 








— 


£- 


-V 




r < e 




= 


i- 


- a 




r > e 



(50) 



There are three types of orbits possible: those which are fully inside the disc, those 
which are fully outside the disc and those which go both inside and outside the disc. In 
the limit of e going to zero, the first type of orbits are clearly irrelevant. It is easy to 
see that for the second type of orbits one must have i ^ a. These are insensitive to the 
flux in the limit and are thus identical to the orbits of the oscillator. 

For the last type, an interior turning point is possible only for < £ < a. The 
condition that such an orbit must also have an exterior turning point limits £ to a. We 
are interested in computing the change in the angular coordinate from the entry into 
the disk till exit from it. Explicit computation shows that the change in the angular 
coordinate goes to zero as e goes to zero. Thus such an orbit reflects at the origin. Note 
that these are precisely the radial orbits. 

To summarize, a regulated classical modeling for anyons is generically stipulated by 
giving the behaviour of the flux near the origin. It amounts to cutting out a disk of 
radius e and filling it up with non-singular fields. The classical non-radial orbits then 
are exactly same as those of the oscillator except for the replacement I ^ I — a. The 
radial orbits reflect at the the origin and are termed half orbits since their period is half 
of that for the other orbits. 
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